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A semiconductor quantum dot can generate highly coherent and indistinguishable single photons. 
However, intrinsic semiconductor dephasing mechanisms can reduce the visibility of two-photon 
interference. For an electron in a quantum dot, a fundamental dephasing process is the hyperfine 
interaction with the nuclear spin bath. Here we directly probe the consequence of the fluctuating 
nuclear spins on the elastic and inelastic scattered photon spectra from a resident electron in a 
single dot. We find the nuclear spin fluctuations lead to detuned Raman scattered photons which 
are distinguishable from both the elastic and incoherent components of the resonance fluorescence. 
This significantly reduces two-photon interference visibility. However, we demonstrate successful 
screening of the nuclear spin noise which enables the generation of coherent single photons that 
exhibit high visibility two-photon interference. 


Indistinguishable single photons are an essential re¬ 
source for quantum photonic logic gates and network¬ 
ing m- Among the various approaches to generate iden¬ 
tical light quanta, resonance fluorescence (RE) from a 
semiconductor quantum dot (QD) [5H5] is one of the 
most promising for practical technological implementa¬ 
tion. The RF spectrum is composed of elastic and in¬ 
elastic scattered light [THin]. The elastic or Rayleigh 
scattered light, first measured in homodyne absorption 
experiments on a QD HU, has the first order coherence 
properties of the laser but the second order coherence 
properties of the emitter [T^. This light, which can be 
imprinted with an arbitrary phase or temporal profile [3] , 
is fundamentally indistinguishable. RF can therefore re¬ 
lax the requirement for a perfectly stable, transform- 
limited optical transition that is difficult to realize in 
the solid-state. In particular, compared to non-resonant 
excitation followed by spontaneous emission, RF helps 
overcome the relatively slow spectral fluctuations caused 
by charge noise in the QD environment [iniin]. 

In addition to the elastic and incoherent components 
of RF, near resonant excitation can lead to Raman scat¬ 
tering into another ground state. Raman scattering from 
a QD with a resident electron or hole can have several 
attractive features. First, below saturation, the coher¬ 
ence of the Raman photons is determined by the ground 
state dephasing rather than that of the excited state [Ti¬ 
lls]. Second, Raman photons are highly tunable: their 
energy is determined by the detuning of the driving 
field rather than the fundamental optical transition en¬ 
ergy [T1|TS1[II]. Finally, due to clean selection rules in 
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a QD, the polarization of the Raman scattered photons 
can be linked to the spin of the hnal state [TSTHSj . en¬ 
abling spin-photon entanglement |19II21] and raising the 
prospect for quantum networks [52|. Typically, spin-flip 
Raman photons are generated using an in-plane external 
magnetic field to mix the spin states [23] . 

A III-V QD also contains an intrinsic effective mag¬ 
netic field, the Overhauser field {Bn), due to coupling 
of the QD’s constituent atoms’ nuclear spins with the 
electron spin via the contact part of the hyperfine inter¬ 
action. Incomplete cancellation of the finite random nu¬ 
clear spin orientations leads to an effective field Bn ^ 10 
- 30 mT that fluctuates on a 10“^ s timescale [Hi^ISS]. 
The electron spin precession around Bn leads to ensem¬ 
ble dephasing on ns timescales [24H26] . The fluctuating 
arbitrary orientation of Bn also leads to an almost al¬ 
ways present in-plane component i?j(f which affects the 
optical properties of a quantum dot. Kuhlmann et al. 
have investigated the role of nuclear spin noise on the 
neutral (Y^) and negatively charged (Y^“) states of a 
QD while Hansom et al. recently exploited to 
achieve coherent control and two-color coherent popula¬ 
tion trapping with X^~ [28] . Here we report new insights 
on the effect of B^j,^ on the RF spectrum of a single QD 
and its consequences on indistinguishable photon gen¬ 
eration. We find spin-mixing causes Raman scattered 
photons clearly distinguishable from the elastic spectrum 
that reduce the visibility {ly) of two-photon interference 
(TPI) for X^~. We demonstrate successful screening of 
B^j,f for both A° and X^~ via either an effective magnetic 
field due to exchange interaction (Jq) or a perpendicu¬ 
lar external magnetic field (S/j.^), respectively. In each 
case, the spectra exhibit near ideal two-level behaviour 
and ultra-high visibility TPI is achieved. 

The self-assembled InGaAs QDs studied in this work 
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FIG. 1. Left: (a) Level diagrams for and X^~ at Be^t = 0 T, and X^~ at > 0 T, where 1+) and |—) are the symmetric 
and anti-symmetric neutral exciton eigenstates, respectively; 5o is the electron-hole exchange interaction energy, and the green 
arrows represent the driving fields, (b) Plot of the ratio of elastic scattered photons to the total intensity. The data points are 
based on the fits to the RF spectra and the curve is a fit using Equation 1 in the main text. Right: high-resolution spectra of 
(c) X°, and (d) X^~ at B^^t = 0 T, and (e) X^~ at B^^^ = 0.6 T, above (If > fls), near (II ~ fig) and below (fl < fig) the 
saturation Rabi frequency. The green (red) line is a single (multiple) Lorentzian fit to the elastic (inelastic) component of the 
RF. The black curve is the total Lorentzian fit. 


were embedded in a GaAs Schottky diode for determin¬ 
istic charge state control. Numerous QDs from two sam¬ 
ples at T = 4K were investigated providing consistent re¬ 
sults. Sample 1, previously used in Ref. [29l yields « 350 
kHz from a single QD on a single photon detector at sat¬ 
uration in CW RF. Sample 2 has the identical Schottky 
diode design but the QDs are positioned at an anti-node 
of a fifth-order planar cavity on top of a Au layer which 
functions both as a mirror and Schottky gate. Sample 2 
yields « 10 times higher count rates due to the planar 
cavity design [3D]. The RF, obtained using orthogonal 
linear polarizers in the excitation and collection arms of 
a confocal microscope I1I1ID1II3I, is characterized by 
three techniques: high-resolution (27 MHz) spectroscopy 
using a Fabry-Perot interferometer (5.5 GHz free spec¬ 
tral range), a Hanbury-Brown and Twiss interferometer 
to measure g^ir), and an unbalanced Mach-Zender (MZ) 
interferometer (Atdeiay = 49.7 ns) with polarization con¬ 
trol in each arm to measure post-selected, two-photon in¬ 
terference nia. In the MZ setup, the beamsplitters have 
nearly perfect 50:50 reflection:transmission. For TPI ex¬ 
periments, a grating (bandwidth f/2 tt ^ 1.9 GHz) is 
used to spectrally filter the zero-phonon line from the 
majority of the phonon sideband. 

Figure reveals our main spectroscopic result. Here, 
we compare the resonantly scattered zero-phonon photon 
spectra from three scenarios: and X^~ with Bext = 0, 


and X^~ with = 0.6 T. The spectra are recorded 
over a range of Rabi frequencies (D), from H <C Dg where 
coherent scattering dominates to H 3> Dg where inco¬ 
herent emission dominates. Here the saturation param¬ 
eter Hg is given by where Ti is the lifetime 

and T 2 the dephasing time. Fig. (c) shows that near 
ideal elastic and inelastic spectra are observed for 
although for D ^ Dg the sidebands of the Mollow triplet 
are slightly broadened from Mollow’s theory due to a 
modest amount of slow spectral fluctuations (which do 
not lead to pure dephasing) m- The central incoher¬ 
ent Lorentzian linewidth is measured as 318 ± 9.3 MHz 
(FWHM), closely matching the Fourier transform limit 
of the measured Ti value of 550 ± 40 ps. Crucially, when 
the population of the excited state is minimal (D ^ Dg), 
the elastic component dominates the spectrum con¬ 
firming minimal pure dephasing. In Fig. (b) we plot 
the ratio of the elastic component to the total spectrum 
(ER). The solid line is a fit given by [7l [3T]: 

where is the ensemble dephasing time obtained in 
time-averaged spectra. By fitting the complete data set 
we can determine the only unknown parameter, : ac¬ 
curately and we find = (1.96 ± 0.08) Ti = 1.08 ± 0.09 
ns for X^. 
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FIG. 2. (a) High resolution spectra as a function of laser 

detuning (A) demonstrate the central and sidebands are de¬ 
pendent on the laser energy, (b) The integrated intensity of 
the detuned spectra reveal a Lorentzian lineshape (red curve), 
(c) A simulated spectrum based on the 4-level model including 
Bn (black curve) fit to the experimental data (black points) 
for Q — 0.25ns. The green line represents elastically scat¬ 
tered photons. The red and blue lines represent the red- and 
blue-detuned Raman photons as shown in the level diagrams 
for two different values, (d): Level diagrams for two dif¬ 
ferent values of For constant excitation energy, changes 
in changes the splitting of the Raman transitions from 
the laser: fluctuations in the Overhauser field therefore result 
in broadened Raman sidebands. 


The spectra for with = 0 in Fig. Bb) 

deviate significantly from the near ideal two-level be¬ 
haviour exhibited by X^. All X^~ spectra show a re¬ 
duced elastically scattered component compared to X°. 
Additionally, a doublet surrounding the elastic peak is 
clearly evident. For the doublet is sepa¬ 

rated by « 190 MHz and each peak in the doublet 
has a linewidth of « 210 MHz. Fig. shows that 
the ER saturates at only « 26 % and the fit yields 
r| = (0.52 ± 0.02) Ti = 0.40 ± 0.07 ns. However, the 
near ideal two-level behaviour can be recovered with a 
modest magnetic field in the growth direction, as shown 
in Fig. [^e) for = 0.6T. Based on the fit to the ER, 


Ti = (1.94 ± 0.01) Ti = 1.51 ± 0.07 ns. 

To investigate the origin of the sidebands for we 
characterize the RF spectrum as a function of laser de¬ 
tuning A for n <C ris. Here we observe that the entire 
spectrum follows A (Fig. [^a)), and the integrated in¬ 
tensity of the spectrum exhibits a Lorentzian lineshape 
with a linewidth of « 850 MHz (Fig. [^b)). This ob¬ 
servation eliminates the possible presence of incoherent 
emission in the spectrum, which is independent of laser 
detuning. We therefore identify the central peak as the 
elastic component and the doublet as spin-flip Raman 
scattering. A schematic of the effect is shown in Fig.j^d), 
in which we depict the four-level system for two values 
of the i?j(r component which split the ground-state elec¬ 
tron spin states. We assume degenerate excited states 
as the hyperhne interaction for hole spins is significantly 
reduced compared to the electron for [35J |331 • The 
solid green arrows represent the driving field on reso¬ 
nance with the unperturbed (e.g. = 0) transitions. 

This driving field determines the elastic scattering spec¬ 
trum. The red and blue dashed lines represent red- and 
blue-Raman scattering detuned by 2/rB<7e7?jl^ = AEz- 
The linewidth of these photons is inversely proportional 
to the incoherent spin-flip time. The experimental mea¬ 
surement time of a few minutes for each spectrum signi¬ 
fies a time-averaged measurement of « 10® orientations 
of Bn, effectively broadening the observed Raman side¬ 
bands according to the variance Notably, we also 

observe QDs which exhibit sidebands with smaller de¬ 
tuning than in Figure 2(a) that can be caused by either 
reduced Bn or gl- Further details on numerical simula¬ 
tions of the resonantly scattered power spectrum based 
on the optical Bloch equations of the 4-level system with 
a linearly polarized driving field are provided in the Sup- 
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FIG. 3. Second-order correlation measurements on the 
Zeeman-split X^~ transitions at = 0.6 T and Q — O.Sfls. 
Single laser excitation of the (a) blue- and (b) red-shifted lines 
show bunching around r = 0 with decay time ~ 20 ns. Simul¬ 
taneous two-laser excitation (c) shows suppressed bunching, 
demonstrating that its origin is spin pumping. 
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FIG. 4. Second-order correlation and two-photon interference measurements on (a) (b) X^~ at Bext = 0 T, and (c) X^~ 

at = 0.6 T. Blue (red) solid lines are (deconvolved) fits to the data. Near perfect TPI visibility is observed for X^ and 
X^~ with screening of nuclear field fluctuations. For X^~ at Be^t = 0 T, i/ is reduced to ~ 0.67 due to Bj^. 


plemental Information. Numerical simulations to fit the 
experimental spectra yield a mean energy of = 

264 ± 113 MHz, with a standard deviation gusSBN of 
120 ± 6 MHz. The black curve in Figure (c) shows 
this simulation. Assuming a typical isotropic electron g- 
factor of -0.6 corresponds to a mean Overhauser field Hjv 
= 31 ± 13 mT and SBn = 14 ± 1 mT. 

As the Overhauser field is ubiquitous in HI-V QDs, 
the spectra from both X^~ with modest and A° 

demonstrate successful screening of the nuclear spin fluc¬ 
tuations. In the case of the exchange interaction 
(i5o) generates a fine-structure splitting of « 10 /reV, 
much larger than gsgeB^L = 0.06 /reV. In the case of 
Bn for X^~, stabilizes the electron spin 

which then precesses around Btot = B^^^ + Bn pS] . 
This significantly reduces the branching ratio for spin-flip 
scattering, which can be estimated by (Bn )"^) 
[34] . Hence, the nuclear spin fluctuations are effectively 
screened and Raman scattered photons are not visible in 
the spectra of Fig. [^e). 

As a precursor to TPI experiments, we first investigate 
the effect of on the second-order coherence (g^(r)). 
Fig. I^a) shows the auto-correlation from A° with a raw 
value of g^{0) = 0.20 ± 0.02, limited by detector jitter. 
Upon deconvolution of the detector response, g^{0) —>■ 0. 
For X^~, while near ideal anti-bunching is again demon¬ 
strated, significant bunching around the r = 0 dip is ob¬ 
served at = 0.6 T (see Fig. |^(a)-(b)) due to optical 
spin pumping caused by electron spin mixing [531HH El] ■ 


The bunching amplitude is determined by the spin pump¬ 
ing fidelity and the decay of the bunching is determined 
by the spin-relaxation dynamics. A two-color excitation 
scheme can be used to frustrate the spin pumping mech¬ 
anism, which unambiguously demonstrates that the na¬ 
ture of the bunching is spin initialization (Fig. ic)). Al¬ 
though we have observed small bunching amplitudes even 
at Bext = 0 for a few QDs, typically it is negligible for 
on resonance driving as in Fig. El(b). 

The canonical test for photon indistinguishability is 
two-photon interference, and the visibility is defined as 
^ = [5i(0) - 5 j|( 0 )]/gi( 0 ), where gl(0) and 5 ^( 0 ) use 
orthogonal and parallel linear polarizations, respectively, 
in the two MZ arms. TPI measurements from A° at 
Q = 0.5Ds (Fig. Qa)) reveal near perfect photon in¬ 
distinguishability, with a deconvolved (raw) fit value of 
ly = 0.99 ± 0.02 = 0.58 ± 0.04). This impressive re¬ 

sult is expected based on the near-ideal high-resolution 
spectrum, confirming that 6o screens Bn for A°. In con¬ 
trast, the same measurement on the X^~ with Bf,xt = 0 
T (Fig. I^b)) shows reduced v: a deconvolved (raw) fit 
value of i/ = 0.67 ± 0.04 {i/ = 0.50 ± 0.05) is measured. 
Nevertheless, a high degree of indistinguishability is still 
observed, somewhat surprising considering the spectra 
in Fig. m- However, in the ’frozen-fluctuation’ model, 
Bn is static within the Atdeiay = 49.7 ns, signifying that 
only one spin-flip Raman transition energy is relevant 
for TPI. In other words, only narrowband Raman pho¬ 
tons will accompany the elastic and modest incoherent 
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components of the RF spectrum over this timescale at 
n = 0.5 fig. Finally, applying = 0.6 T to screen the 
ground state electron from the nuclear field fluctuations 
recovers high visibility TPI with a deconvolved (raw) fit 
value of z/ = 0.97 ± 0.03 {v = 0.64 ± 0.06)(see Fig. |^c)). 
This near perfect indistinguishability meets the expecta¬ 
tion based on the spectra in Fig. [^e). 

In summary, we have characterized the effect of nu¬ 
clear field fluctuations on the RF spectrum and photon 
indistinguishability from and X^~ in self-assembled 
InGaAs QDs. For the A°, the large Sq relative to Bj^ 
suppresses the effect of nuclear field fluctuations even at 
Bext = 0 T. For a QD with a resident electron in the 
ground state, the RF spectrum exhibits an inelastic dou¬ 


blet due to spin-flip Raman scattering which reduces the 
TPI visibility. Application of B^^^ screens the electron 
spin from B^^, drastically decreasing the rate of Raman 
scattering and recovering near-ideal indistinguishability. 
An interesting prospect would be to investigate the pos¬ 
itively charged exciton X^'^, for which the ground hole 
spin states can be more robust against nuclear spin noise 
at B^xt — 0- 
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Supplementary Information: Screening nuclear field fiuctuations in quantum dots for 

indistinguishable photon generation 


In this supporting document we provide the following: details on the determination of two- 
photon interference visibility, modelling of the bunching due observed in second-order corre¬ 
lations due to optical spin pumping of the charged exciton, and a theoretical model for the 
resonance fluorescence spectrum of the transition under non-zero nuclear field. This model 
provides the fitted spectrum shown in Fig. 2(c) of the main text. 


SI. DETERMINATION OF TWO-PHOTON 
INTERFERENCE VISIBILITY 

A. Deconvolution 


Determining the indistinguishability based on the ex¬ 
perimental data is achieved using the canonical functions 
for two-photon interference (TPI) [IHS], assuming 50% 
reflection and transmission at both beamsplitters: 

5'(T) = l-e-l"l/"i (SI) 

A'r)+5'^(T + Ar))-b^5^(T) (S2) 

= ^(g^(r-Ar)-b5f2(T-_gAr))(l-i;e"2|^l/^‘=) 

+ (S3) 


where ti is the exciton fluorescence lifetime, At is the 
fibre delay in the Mach-Zehnder interferometer, v is the 
visibility and Tc is the photon coherence time. 

These functions are then convolved with the set-up in¬ 
strument response function (IRF), which is determined to 
be well approximated by a Gaussian peak with full width 
at half maximum (FWHM) of 750 ps. This yields an an¬ 
alytical expression for the convolved function in terms of 
real physical parameters. We use this convolved func¬ 
tion to fit the raw data, and the deconvolved fit can be 
plotted by using these parameters extracted from the fit 
in Eqs. ( Slf^ ). Ti is extracted from a fit to the second 
order coherence function using ^^(t), and is then used as 
a fixed parameter in the TPI fitting. 

In the convolved TPI data, there is a certain amount of 
ambiguity between the effects of visibility and coherence 
time on raw dip depth that determines the photon indis¬ 
tinguishability. As photons become more distinguishable, 
V drops from 1 to 0 for perfectly distinguishable photons 
(e.g. photons with perpendicular polarisations, photons 
with very different energies), which reduces the depth of 
the central dip for (r) even after deconvolution. How¬ 
ever, when the interacting photons are indistinguishable 
but have very short coherence times, the width of the 
central dip for (r) from 50% to 0 will significantly nar¬ 


row. In this case the experiment is more strongly affected 
by the system IRF [3]. This raises the question whether 
a dip depth above zero in the raw data is an effect of 


photon distinguishability or very short coherence time. 
Fig. [ST] demonstrates this ambiguity. 

In this manuscript, near ideal transform limited 
linewidths are observed in the high resolution spectra. In 
the case of X^~, we show the non-ideal behaviour is due 
to Raman scattering at energies detuned from the cen¬ 
tral transform limited peak in the spectra. We therefore 
assume this is the origin of the reduced raw (experimen¬ 
tal) TPI visibility, rather than reduced coherence or pure 
dephasing. 


B. Spin pumping 


As observed in Fig. 3, spin pumping in the 4-level sys¬ 
tem leads to bunching in ^^(t) and TPI plots. This is 
captured in the fitting process by r epla cing t he g ^jr) ex¬ 
pression given by Eq. (SI) in Eqs. (S2| and (p3| with an 


exponential dip multiplied by an exponential peak with 
background 1 and a longer decay time t;,: 


gl{T) = (1 - e-l"l/^i)(l + fee-'"'/"'’) 


(S4) 


where b is the amplitude of the bunching. In Fig. 3, this 
yields b = 0.85 ± 0.01 and Tf, = 15.9 ± 0.3 ns. 


SII. RESONANCE FLUORESCENCE 
SPECTRUM OF THE TRANSITION UNDER 
NON-ZERO NUCLEAR FIELD 

The electron and hole wavefunctions of self-assembled 
III-V quantum dots (QD) extend over a large number of 
spin bearing nuclei. The nuclear spins generate an ef¬ 
fective magnetic field, called the Overhauser field, that 
couple with the spins via the hyperfine interaction and 
can lead to relaxation and dephasing. In the context 
of the low temperature experiments performed in this 
manuscript, the “frozen-fluctuation model” |3], in which 
the nuclear spins are considered to be stationary for a 
timescale (« 1 fis) much greater than the optical transi¬ 
tion lifetime (« 1 ns), is most relevant for self-assembled 
quantum dots In this picture, the electron spin pre- 

cesses around the hyperfine field on a ns timescale. Here 
we are interested in changes to the energy levels and se¬ 
lection rules caused by the random Overhauser field ori¬ 
entation. In particular, we develop a theoretical model 
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FIG. SI. Simulated TPI and visibility plots for two different situations showing similar convolved shapes. Red: plot given by 
Eqn. S3 Blue: plot given by convolving Eqn. S3 with a Gaussian IRE with FWHM 700 ps. n = 700 ps, At — 15 ns. (a): TPI 
data with good coherence (tc = 2ti) but poor visibility (v = 0.6). (b): TPI data with perfect indistinguishability (v = 1) but 
poor coherence (tc = 0.5 ri). (c) & (d): Plots of visibility for the cases in (a) and (b), respectively. 


for predicting the resonance fluorescence spectrum of the 
photons scattered from the transition of the QD. 
Our model features four levels, two of those represent¬ 
ing the electron spin of the QD’s ground state, and the 
other two the unpaired (heavy) hole spin states of the 
first excited state. 

We let the electron spin levels be Zeeman split by the 
total magnetic field, i.e. the sum of the externally applied 
field Bg^t plus the effective Overhauser contribution Bj^i 
(in the frozen fluctuation regime). By contrast, we as¬ 
sume the hole spin only ‘sees’ the external field that is 
applied along the z-direction. (Due to the reduced hole 
Bloch functions at the nuclear positions, the hyperfine in¬ 
teraction coupling coefficient between a hole and nuclear 
spins is only « 10 % that of the electron spin coupling 
coefficient BM- Further, the hyperfine interaction for 
an ideal heavy hole takes on an Ising-like form such that 
heavy hole spin dephasing is greatly suppressed by an 
in-plane magnetic field [10]). As the basis for our model 
we adopt the Zeeman eigenstates of both electron and 
hole (but with respect to different magnetic fields as dis¬ 
cussed above), see Fig. S2 for a schematic depiction of 


the relevant level structure. 

Writing the total magnetic field B = Bg^tZ + i?jv as 

B = B [sin 6 cos (j)x + sin 6 sin (j)y + cos 9z] (S5) 

the electron spin Zeeman [pi ■ B) eigenstates are 


1+) = cos (^0 It) + sin ^0 |t) , (S6) 

I-) = sin It) - 6*“^ cos It) . (S7) 


Denoting the trion states with hole spin (anti)parallel 
to the z-direction as (|1|)) |11'),[II] we obtain the following 



















optical dipole transition matrix elements: 



FIG. S2. The X^~ transition under an Overhauser field with 
an in-plane component. When the total magnetic field is not 
completely parallel to the growth (z) direction the electron 
spin-mixing alters the selection rules and all four transitions 
become dipole allowed. 


(tr|r |-h) oc COS c_ , (S8) 

(JJ.|r |-h) oc e*'^sin c+ , (S9) 

(tr|r |-) oc sin c_ , (SIO) 

(JJ-I ^|-) cx cos c+ , (Sll) 

where c± refers to the two possible circular polarisation 
states. Including the dependence of the selection rules on 
the magnetic field direction and under a rotating wave 
approximation, the system Hamiltonian is then given by 


^ l+)(+l +1 l-)(-l + (^- y) + (^ + y) 

+ (e*‘^e-c-h|+)(l||+i?.c.)-h^cos^0 c_ |-h)(tr| + i?.c.) 

- ^ cos (e*'^e • c+ |-)(.ll| -h H.c.) + ^ sin (e • c_ |-)(iy| -h H.c.) , (S12) 


where H.c. denotes the Hermitian conjugate, 5^ = 
gelJ-sB, Sh = ghfJ-sBext, is the optical Rabi frequency 
and A the laser detuning from the QD transition in the 
absence of any magnetic field, and e is a unit vector in 
the direction of the polarization of the driving laser field. 


The selection rules Eqs. ( S^Sll ) also imply sponta¬ 
neous emission decay channels with polar angle 9 depen¬ 
dence as follows 


r+^ = rcos2Q) , (S13) 

r+4 = rsin2Q^ , (S14) 

T-f = Tsin^(J^ , (S15) 

T_^ = Tcos^(J^ , (S16) 


where T = l/Ti is the inverse lifetime of the QD. These 
processes are depicted schematically in Fig. |S2[ Further, 
we define a lowering operator for each circular polarisa¬ 


tion 

ai^^ = (x/r;il+)(^| + y/rli|-)(^|)c+, (S17) 
ai^) = (v^|+)(^| + V^|-)(^l)c- . (S18) 

Using these equations, we can obtain the temporal evo¬ 
lution of the density matrix, given by numerically inte¬ 
grating the master equation 

^ p] + L{a^}^)p + L{a^}^)p , (S19) 

where L{o)p = op ol — (ol o p -|- p o'! o)/2 is the standard 
Lindbladian dissiDator|12). 

We proceed to calculate the first-order correlation 
function = (cr+(t-|- T)cr_(t)), which will give us 

access to the fluorescence spectrum of the photons scat¬ 
tered by this four-level system. In our case we have (due 
to the two orthogonal polarisations) 

G«(r) = + + (af + T)aL^^(t)) . 

(S20) 

The fluorescence spectrum is independent of both the 
azimuthal angle (j) and the polarisation angle for a linearly 
polarised driving field. It is calculated by taking the real 
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part of the Fourier transform of the first-order correlation 
function. For all results presented below, the 5-peaks of 
elastically scattered photons and the Raman transitions 
are artificially broadened by making the expected number 
of photons equal to the area of a Lorentzian curve of 
width 27.5 MHz, corresponding to the resolution of the 
Fabry-Perot interferometer used in the experiments. 

The (elastically) scattered Rayleigh peak is present for 
any orientation of the magnetic field, but its intensity 
depends on 9. In Faraday geometry {6 = 0), there is 
no Raman transition and the elastic peak is accompa¬ 
nied by two detuned side peaks at ±5e/2 which possess 
a width limited by the decay rate of the optical transi¬ 
tion [Fig. S3 (a)]. The relative amplitudes of the central 
to detuned peaks is determined by the Rabi frequency. 
When the magnetic field is between the Faraday and the 
Voigt geometry, the relative weight of the two vertical 
transitions decreases whilst the two diagonal transitions 
become allowed. This gives rise to two Raman transi¬ 
tions [Fig. [^(b)]. These appear as 5-peaks, but have 
here been broadened by the detector resolution. In the 
Voigt geometry (0 = 7r/2), all four transitions have equal 


weight, thus maximising the intensity of the Raman tran¬ 
sitions, which now account for 50% of the total number 


of scattered photons [Fig. S3 (c)]. 

To reproduce the long time measurements, it is nec¬ 
essary to average several spectra with different instan¬ 
taneous nuclear fields Bj^. Considering that the Over- 
hauser field has no preferential direction, we will realise 
(Hat) = 0 (as required) upon performing uniform aver¬ 
aging over the solid angle. It remains to specify a radial 
distribution of to represent the magnitudes of the in¬ 
stantaneous Overhauser fields. We consider the case of a 
normal distribution with generally finite mean (Hat) and 
variance 6Bn. For (Bn) = 0 the difference between the 
energy of the Raman scattered and the Rayleigh photons 
is typically small, and consequently the average over the 
Raman peaks yields a single Gaussian envelope of width 
proportional to SBn [Fig. [^(a)]. By contrast a finite 
mean (B^) >0 separates the Rayleigh from the Raman 
peaks, leading to two shifted Gaussian Raman envelopes 
either side of the central elastic peak [Fig. S4 (b)]. This 
latter case gives very close agreement with our measured 
spectra shown in the main text [see also Fig. |S4| (c)]. 
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frequency (GHz) 


FIG. S3. Resonance fluorescence spectra with flxed Overhauser fields for B^xt = 0 MHz: (a) Faraday geometry with Rabi 
frequency H = O.lHsat. At sufficiently strong driving (here ten times stronger than in the other two panels of this figure), 
the central elastic peak is accompanied by two inelastic peaks detuned by ±5e/2 {Bn = 500 MHz), originating from radiative 
decay with orthogonal polarisation to that of the excitation photon, (b) 6 = 7r/4 and H = O.OlHsat. The two side peaks at ±(5e 
{Bn = 250 MHz) are due to Raman scattered photons (broadened to the interferometer resolution), (c) Voigt geometry with 
other parameters as in (b): the Raman transitions reach maximum intensity in this configuration. 



frequency (GHz) 


FIG. S4. Resonance fluorescence spectrum averaged over a fluctuating Overhauser field, using F = 1.33 GHz, H = O.OlHaat, 
Bsxt = 0 MHz, and an isotropic angular distribution of the Overhauser field throughout, (a) {Bn) = 0 MHz and SBn = 200 
MHz: the Raman peak envelope forms a single Gaussian with width proportional to SBn- (b) {Bn) = 300 MHz and 5Bn = 100 
MHz: the Raman peak envelopes can be resolved separately either side of the central elastic peak, (c) Our averaged full 
simulation of the spectrum with {Bn) = 264 MHz and 5Bn = 113 MHz provides an equally good fit to the experimental data 
as directly fitting (i.e. without underlying calculation of the spectrum) two Gaussian peaks with mean 0.208 ± 0.005 MHz and 
standard deviation 0.114 ± 0.006 MHz plus a Lorentzian peak with mean equal to 0 MHz and width equal to 27.5 MHz (FPI 
resolution). 







































